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Reducing the impact of loss-of-control conditions on commercial transport aircraft is a 
primary goal of the NASA Aviation Safety Program. One aspect in developing the 
supporting technologies is to improve the aerodynamic models that represent these adverse 
conditions. Aerodynamic models appropriate for loss of control conditions require a more 
general mathematical representation to predict nonlinear unsteady behaviors. In this paper, 
a more general mathematical model is proposed for the subscale NASA Generic Transport 
Model (GTM) that covers both low and high angles of attack. Particular attention is devoted 
to the stall region where full-scale transports have demonstrated a tendency for roll 
instability. The complete aerodynamic model was estimated from dynamic wind-tunnel data. 
Advanced computational methods are used to improve understanding and visualize the flow 
physics within the region where roll instability is a factor. 
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Nomenclature 


Fourier coefficients 

deficiency function parameters 

wing span, ft 

rolling -moment 

yawing-moment coefficient 

side-force coefficient 

mean aerodynamic chord, ft 

deficiency functions 

frequency, Hz 

reduced frequency, 7rbf / V 

no. of harmonics in Fourier expansion 

number of data points 

roll and yaw rates, rad/sec 

coefficient of determination 

reference area, ft 2 

estimated standard error 
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dimensional time constant, sec 
time, sec 
velocity, fps 

angle of attack, rad or deg 

initial angle of attack in forced- 

oscillation experiments, rad or deg 

sideslip angle, rad or deg 

state variable 

standard error 

dummy integration variable 


non-dimensional time constant, 



roll angle, rad or deg 
yaw angle, rad or deg 
angular frequency, rad/sec 
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Subscripts 

A = amplitude 

a = l, n, or Y 

E = measured value 

Superscripts 

A = estimated value 

= mean value 


Aerodynamic Derivatives 

C n - = — — where a = l,n, or Y 

a P dp 

= b_dC< L b dCg 

a P 2V dp ’ “ r 2V dr 


I. Introduction 

R educing the impact of loss-of-control conditions on commercial transport aircraft is a primary goal of the 
NASA Aviation Safety Program. In support of this goal, both analytical and experimental technologies 
continue to be developed to improve the aerodynamic models that represent these adverse conditions. Aerodynamic 
models appropriate for loss-of-control conditions require a more general mathematical representation to predict 
nonlinear unsteady behaviors. Researchers are addressing this problem by developing analytical tools in 
Computational Fluid Dynamics (CFD) and experimental tools using direct measurement of aircraft responses in 
flight and in wind tunnels along with System Identification (SID) technology to determine mathematical models. 

A number of experimental and analytical studies have been performed to obtain aerodynamic models for the 
subscale NASA Generic Transport Model (GTM) and this work has produced significant improvements in 
simulating transport aircraft responses in loss-of-control conditions. Early studies to improve GTM simulation in 
upset conditions were reported by Cunningham et al. 1 in 2004 and Foster et al. 2 * in 2005. In these studies an 
expanded aerodynamic data base was produced and both static and dynamic roll-axis instabilities and nonlinearities 
were identified at stall angles of attack. Later in 2010, Belcastro and Foster 1 provided a study showing that loss of 
control is still the leading cause of accidents worldwide. In 2011, Advani and Field 4 highlighted the need for 
integrated upset prevention and recovery training in flight simulators. 

With this motivation, both computational and experimental research in modeling transport aircraft has continued. 
A recent paper by Morelli 5 in 2012 produced global aerodynamic models for GTM using flight test data. The results 
were obtained using automated perturbation inputs and were presented over a large range of angle of attack. 
Instability was found and characterized as a roll-off departure at 14° angle of attack. These results are consistent 
with the modeling work by Murphy and Klein 6 where dynamic roll instability was found near 12° angle of attack. In 
Ref 6 results were obtained using advanced modeling methods applied to dynamic forced-oscillation tests of the 
GTM. In this later study, nonlinear unsteady responses were identified in the stall regime. The strength of this roll 
instability and the abrupt nature of the roll-off departure warranted further investigation. 

In this study both CFD and SID approaches were used to investigate the subscale NASA GTM roll-axis 
aerodynamics. The aerodynamic model was estimated from dynamic wind-tunnel data and additional insights into 
the flow physics were obtained using advanced computational methods where roll instability and unsteady behaviors 
are a factor. The SID approach is discussed by presenting a simplified version of a more general formulation of the 
transport aerodynamic model structure using indicial models 7 . The identification problem is pursued by first 
performing harmonic analysis to locate and assess the severity of nonlinear and unsteady behaviors. This analysis 
helped identify appropriate model structure for parameter estimation. Although the global model for rolling moment 
is nonlinear, local models identified by dynamic tests must be evaluated to determine their adequacy as linear or 
nonlinear representations over the range tested. In this study only local models are developed with appropriate 
model structures to reflect the degree of nonlinearity and unsteadiness. These models, however, can form a 
foundation for a global model in simulation. Final parameter estimation was accomplished by using output-error, as 
described in Ref. 8. Data filtering and its impact on harmonic analysis are discussed in detail to clarify key modeling 
steps that lead to determination of linear or nonlinear model structures. Identifying unsteady aerodynamic models 
can be performed efficiently with advanced testing methods that use wide-band inputs 9,10 ; however, for this study 
only conventional single-axis, single-frequency, forced-oscillation data are available from previous wind-tunnel 
studies. 

Advances in computational methods offer an opportunity to complement the wind-tunnel data with physics- 
based analyses that effectively provide an “analytical” wind tunnel. These methods allow the use of advanced test 
techniques and motions not readily available when applying conventional experimental methods and thus allow 
more flexibility in understanding and estimating key aerodynamic modeling parameters. The combined experimental 
and analytical experiments offer an opportunity to gain greater insights into aerodynamic phenomena observed in 
flight. An initial investigation demonstrating the potential benefits of this approach was reported in Ref. 1 1 , where 
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linear unsteady aerodynamic models were identified from both wind-tunnel measurements and CFD simulations. 
Aerodynamic model comparisons were obtained using the NATO RTO AVT-161 SACCON 12 wind-tunnel model. 

II. Wind Tunnel Measured Data 

Wind-tunnel tests were conducted at the NASA Langley 14x22 Wind Tunnel in 2009 to obtain lateral-directional 
data appropriate for identification of unsteady aerodynamic models. Analysis of this data was reported in Ref. 6, 
where analysis of the yaw-axis measurements was the primary focus. A 5.5% scale model, representing a 
conventional twin-engine commercial transport, was tested. Figure 1 shows the basic geometry of the model. Roll 
and yaw forced oscillation experiments were performed over a wide range of frequencies and amplitudes to facilitate 
unsteady model identification. For the roll case, the tunnel speed was 92 fps and measurements were made for roll 
oscillations at angles of attack from -5 to 75 degrees; five frequencies of 0.06, 0.12, 0.23, 0.46, and 0.92 Hz; and 
four amplitudes of 5, 10, 20, and 30 degrees. The frequencies chosen correspond to the non-dimensional 
frequencies: k=[0.015 0.028 0.054 0.108 0.215]. During the test runs, data were obtained from 7 oscillation cycles 
for low frequency data and up to 44 cycles for the remaining data. The sample rate was 250 Hz with a low-pass 100 
Hz analog filter. The resulting data were further filtered with a 4 Hz low-pass digital filter to remove unwanted 
frequency content. The 4 Hz filter was run in both directions to ensure no phase error was added to the data. 

III. CFD Simulated Data 

CFD predictions were performed with the USM3D flow solver 13 that is part of the NASA Tetrahedral 
Unstructured Software System (TetrUSS) 14 . USM3D is a parallelized tetrahedral cell-centered, finite volume 
compressible Reynolds-averaged Navier-Stokes (RANS) flow solver. The term “cell centered” means that the finite 
volume flow solution is solved at the centroid of each tetrahedral cell. Inviscid flux quantities are computed across 
each tetrahedral cell face using various upwind schemes. Spatial discretization is accomplished by a novel 
reconstruction process, based on an analytical formulation for computing solution gradients within tetrahedral cells. 
The solution can be advanced in time by a second order “physical” time step scheme, a second order “dual” time 
step scheme, or to a steady-state condition by an implicit backward-Euler scheme. Several turbulence models are 
available including the one-equation Spalart-Allmaras (SA) model and several two-equation models. The two- 
equation models available are the Jones and Launder k-s model, Menter Shear Stress Transport (SST) model, 
nonlinear Algebraic Reynolds Stress Models (ARSM) of Girimaji and Shih/Zhu/Lumley, and the Wilcox 1988 k-co 
model. Detached Eddy Simulation (DES) has been implemented in all of the turbulence models. A capability to trip 
the flow at specified locations on aerodynamic surfaces has been implemented for the k-s turbulence model, but 
fully turbulent flow was assumed for the results in this paper. USM3D has capabilities for overset grids and dynamic 
grid motion, the latter being used in the current study. Published guidelines 15 for computing dynamic forced 
oscillation solutions were employed for the CFD computations in this study. CFD simulations were generated that 
correspond to the roll forced-oscillation wind-tunnel tests performed on the GTM in the regions of interest. 

IV. Aerodynamic Modeling 

Simplified model equations for the lateral coefficients Q , C n , and C Y represented by C a , where a = l, n, or Y, 
can be developed from a general form of the indicial model equations 7,16 ' 18 . Each coefficient is considered in the 
form 

C a (0 = C a (0) + I C a B tt “ r)/3(z)dz + 2- J C a (t ~ z)p(z)dz + ^-\C (t- z)r(z)dz (1) 

0 H 2v o y 2y 0 

where C a ^(t), C a ^(t ) , and C ar (t ) are the indicial functions and C a ( 0) is the initial value of C a . Two 

assumptions were adopted to simplify the model: the effect of angular accelerations p and r on any coefficient can 
be neglected and the indicial functions in Eq. (1) can be expressed as 

C a /t) = C a /co)-F a/} (t) ( 2 ) 
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where F a > ( t ) is the deficiency function and C a ^ (oo) is the partial derivative of C a with respect to /? evaluated in 
steady flow conditions. 

The simplified model, which takes into account changes with respect to steady state, has the form 

c a (t) = C ap (<o)fi(t) + A C ap (co) p(t) + A C ar (<o)r(t) - j F ap ( t - x)p(x)dx (3) 

where C a (0) has been subtracted from both sides of Eq. (3). To obtain a model appropriate for identification and 
with a limited number of parameters, the deficiency function is assumed to be a simple exponential function 16 

F ap =ae~'V. (4) 


Models appropriate for an aircraft undergoing one degree of freedom forced oscillation in roll or yaw can be 
obtained using Eqs. (3) and (4). Considering one degree of freedom rolling motion in the tunnel 

C a (t) = C a W\P<t)\ (5) 

where roll angle is related to the sideslip angle by the equation 

Pit) = sin -1 (sin a sin (6) 

Combining Eqs. (3-6), the aerodynamic models can be formulated as 

c a (t) = C (*>)fi(t) + A C (00). p(t) - af (7) 

H F 0 


By introducing 

77 (f) = J P{r)dr 

0 

and applying the Leibnitz integral rule, the state space form of Eq. (7) can be written as 

T)it) = -b l Tjit) + Pit) 


( 8 ) 


(9) 


c a (t) = c ap (co)P(t) + A c ap (CO) p(t) - arj(t). (10) 

From Eq. (7), a steady response can be obtained 16 as 

c a (0 = C ap <t> A sin(or) + C ap </> A k cos (cot) (11) 


where <j> A is the amplitude of roll oscillation, k is reduced frequency, and C ap and C Up are the in-phase and out-of- 
phase components, respectively. These components are related to the model parameters by the equations 17 


C ap =C a p( c °) sin «0' 


X?k 2 

+ rfk 2 


sin a () 


( 12 ) 
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(13) 


C a„ = C a„( co )~ a - 


1 + T[k 2 


sm« n 


where r x is the non-dimensional inverse of b x . 


V. Results and Discussion 

Before identifying the unsteady models described in the last section, analysis is required to determine where, 
with respect to angle of attack and sideslip, the more general models are needed, if at all, and to determine if a 
nonlinear formulation is required. Harmonic analysis along with the coefficient of determination, R 2 , is used for that 
step. Since R 2 is sensitive to noise, the analysis can be made more difficult with low signal/noise data. Before 
presenting model parameter estimates, some discussion is provided to help establish a best practice for harmonic 
analysis and examples are given to demonstrate the effects of noise. Non-zero ordinate values on all plots were 
removed in order to maintain proprietary agreements. 

A. Harmonic Analysis 

A method of harmonic analysis 19 was applied to measured aerodynamic coefficients. A mathematical model for 
these coefficients is 


m m 

C a (t) = Aq+ Y Aj cos [j(9t)+ Y Bj sm( j(ot) a = l,n or Y (14) 

./=> ./=> 

where A {) , Aj , and Bj are the Fourier coefficients. The analysis provides estimates of these coefficients, their 
standard errors, and the coefficient of determination, R 2 . For the model with linear aerodynamics and Aq = 0, the 
aerodynamic in-phase and out-of-phase components can be expressed in terms of the coefficients A 1 and B x . For the 
roll oscillation case the expressions are 


where /? is related to ^by Eq. (6). 

The coefficient of determination, R 2 , indicates the fraction of variation in the measured data explained by the 
variation in the model and is defined as 


C =A 


c„ 


k <t>A 


(15) 


(16) 


R 2 =l-SS E /SS r , 0<R 2 <1 


(17) 


where 


ss E = uc (i)-c a (i )] 2 

i = 1 


(18) 


is the residual sum of squares and 


ss r = z [C a (i)-C a (i)] 2 

i=\ ^ 


(19) 
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is the total sum of squares. C aE (/), C a (/), and C a (/) are the measured, estimated, and mean values, respectively. 

Results of harmonic analysis performed on roll oscillatory data with amplitudes of 10°, 20°, and 30° are 
presented in Fig. 2a, b, and c, respectively. These plots show in-phase and out-of-phase components against angle 
of attack at different frequencies and are plotted with the same scales. The plots indicate very limited frequency 
dependence (unsteady behavior) for the in-phase component for all three amplitudes. Frequency dependence in the 
damping term (out-of-phase component, C tp ) is present at a 0 = 12° for all three amplitudes. For the small 

amplitude case frequency dependence is also present at all a 0 above 24°. For the larger amplitudes the unsteady 
behavior splits into two regions: 24° < a 0 < 30° for amplitude 20°, and 26° < a 0 < 30° for amplitude 30°. Both larger 
amplitudes present some unsteady behavior for a 0 > 40°. Roll damping instabilities are present at a 0 near 12°, 14°, 
and above 30° for the smaller amplitude and at a 0 near 12°, 14°, and above 40° for the two larger amplitudes. As an 
example of yaw, and side force derivatives, similar analyses are shown in Fig. 3 for one amplitude, (j) A = 20°. 
Unsteady behavior and instability occur in approximately the same a 0 locations. 

B. R 2 Analysis 

As follows from Eqs. (17) to (19), the estimates of R 2 are influenced by the value and number of Fourier 
coefficients (harmonic order) in Eq. (14) and the measurement noise in C a (/) . Because the Fourier coefficients are 

mutually orthogonal, the estimates of Aj and Bj will not change with the number of coefficients included in Eq. 

(14). Changes will only appear in the corresponding standard errors and residuals defined by Eq. (18). This makes 
R 2 an effective diagnostic tool to discern the adequacy of a linear first-order model against nonlinear higher 
harmonic models. 

In Fig. 4a the harmonic analysis for the roll moment coefficient is repeated with R 2 shown for the middle 
amplitude case, $ A = 20°. Analysis was performed on data pre-filtered at 4 Hz. At first look, R 2 indicates that 
substantial nonlinarity exists in the stalled flight regions for 12° < a 0 < 40°. However, before making any 
observations on the details of this case, the effects of filtering should be discussed. 

In Fig. 4b the rolling moment coefficient in steady oscillation is presented as mean cycles of measured forced- 
oscillation data and a first order harmonic model. Before computing the mean cycle, three different band-pass filters 
were applied to produce data with substantially different noise levels. The corresponding harmonic models for each 
case are plotted (with same scales) to show the relationship of the data to its estimated linear model. Since a linear 
model can only produce a regular ellipse in this type of plot, 17 any deviations from that regular shape indicate some 
degree of nonlinearity. The plots show that as filtering is done more aggressively, higher-order harmonics are 
removed. Consequently, engineering judgement is needed to avoid filtering too aggressively, removing information 
content and changing the resulting derivative estimates. In this example, the linear model provides a reasonable fit 
for even the least filtered 10 Hz data. 

To further demonstrate the effect of filtering, the best use, and interpretation of R 2 , Figs. 5 and 6 are presented. 
These figures show wind-tunnel measurement data filtered at 10 Hz and 2 Hz for side force and roll moment, 
respectively. In Fig. 5a, R 2 highlights poor adequacy of the linear model primarily for 10° < a 0 < 30°, however, in 
this case filtering at 10 Hz has guaranteed the measurements were dominated by noise (or other responses outside 
the bandwidth of interest, such as structural vibrations) and in turn corrupted R 2 . By filtering down to 2 Hz, as 
shown in Fig. 5b with the same scales, a significant portion of the noise was removed and the R 2 metric now implies 
a linear model is completely adequate over the full range of a 0 . Only at a 0 = 0° are R 2 values still low. This is caused 
by very small side force values at this condition resulting in very poor signal/noise levels. Also, in this case the 
aggressive filtering did not change the in-phase and out-of-phase coefficient values in any significant way. 

The same treatment was performed on roll moment and the results are shown in Fig. 6. In Fig. 6a, with limited 
10 Hz filtering, R 2 implies poor effectiveness of the linear model for 12° < a 0 < 40°. In Fig. 6b, with 2 Hz filtering, 
the low R 2 range has been reduced to 14° < a 0 < 35° and the coefficients have not changed values compared with 
the 10 Hz filtering case. However, in the 2 Hz filtering case a number of conditions are still showing very low R 2 
values. The three lowest R 2 cases at a 0 = 14°, 22°, and 30° are explained by considering measurement and harmonic 
model time histories plotted as response variable versus input variable. These three cases are representative of all the 
low R 2 cases so only these three will be discussed. The first case, at a 0 = 14°, demonstrates steady dynamics but 
with a nonlinear response; the second case, at a 0 = 22°, shows steady dynamics with a linear response; and the third 
case, at a 0 = 30°, presents unsteady dynamics with a linear response. For the set of amplitudes and frequencies 
presented in this study, the fourth possible case, a combination of nonlinear and unsteady responses, was not 
observed. An example of the fourth case was presented in Ref. 6, where a third order harmonic was required to 
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model the cubic response in yawing moment to oscillations in sideslip. The source of the nonlinearity was found to 
be in the static term because of a strong cubic anti -symmetric relationship of yawing moment to sideslip. This type 
of model is common in the lateral-directional axis. In general, further investigation during the parameter estimation 
process is required to determine whether the nonlinearity is required in one or all three major model terms of Eq. (3). 
The three major terms are: (1) static terms, (2) steady-flow dynamic terms, and (3) unsteady terms. 

Figure 7 shows rolling moment and side-force coefficients in steady oscillation for the three a 0 cases with low R 2 
at high and low frequencies and for both side force and roll moment. In order to highlight ellipse shapes the scales 
are not the same. The left side presents plots for the low frequency cases at f - 0.06 Hz, and the right side shows 
high frequency cases at /= 0.92 Hz. At the top of the figure are results at a 0 = 14°, where a low R 2 value occurred 
for the high frequency oscillations. In this case the low R 2 value correctly indicated a nonlinear response since the 
high frequency plot for C/ is an irregular ellipse, and linear aerodynamics can only produce regular ellipses. 
Nonlinear responses are often expected in the lateral -directional axes where a strong cubic nonlinearity typically 
exists with respect to sideslip. This case is also a good example of where engineering judgement is required. 
Although the response at high frequency is nonlinear, respect for the very small values of the derivatives may lead to 
a conclusion that the added complexity of a nonlinear model is not warranted. R 2 does not account for overall 
magnitude of the response variables. The linear responses predicted by R 2 for C Y at both high and low frequencies 
are also confirmed. 

For the second case at a 0 = 22°, the low R 2 value occurred at the lowest frequency. In this case noise has reduced 
R 2 . Plots in the middle of Fig. 7 apply to this case and show that linear models match the measured data reasonably 
well except for the low frequency case of the roll moment, where the response measurement has substantial noise in 
spite of aggressive 2 Hz filtering. In this case, although the noise has reduced the R 2 metric, the estimated linear 
model is still adequate. 

The third case at a 0 = 30°, presents a case where an unsteady model is required and a low R 2 value occurred at 
the lowest frequency (as shown in Fig. 6). Plots at the bottom of Fig. 7 apply to this case and the plots show linear 
models match the measured data reasonably well except for the low frequency case of the roll moment where the 
response measurement has substantial noise. Once again, although the noise has corrupted the R 2 metric, the 
estimated linear model is still adequate. Although the lag in response is not obvious, comparing the second and third 
cases, Fig. 6 highlights the unsteady model requirement due to frequency dependence for the third case at a 0 = 30°. 

Figures 5, 6, and 7 highlight three cautions when using R 2 : (1) the metric R 2 is sensitive to signal-to-noise levels; 
(2) the metric R 2 does not account for the magnitude of response variables, so engineering judgement is required to 
evaluate if the response levels should be modeled; (3) model derivatives should not change values with appropriate 
filtering; consequently sensitivity of derivative estimates should be checked against the degree of digital filtering. 
Maintaining a factor of 5 between the digital filter cutoff frequency and the unsteady model pole locations has 
proven useful. Although filtering at 2 Hz did not change the in-phase and out-of-phase coefficients for this study, in 
practice 4 Hz has been found to be a good initial choice for filtering forced-oscillation data of rigid-body transport 
aircraft. 

C. CFD Simulation Analysis 

A previous study 11 presented initial efforts to compare unsteady model parameters estimated using both CFD and 
dynamic wind-tunnel data. The promising results of that study have prompted the authors to extend the use of CFD 
simulation in this study. Although it is not the objective of this study to validate CFD, a number of additional 
insights into the physics behind the wind-tunnel results can be inferred, given the reasonably comparable 
aerodynamic predictions between CFD and wind-tunnel data. Further work is required to validate CFD applied to 
dynamic testing; however, with the use of CFD best practices, 15 researchers may obtain corroborating data or 
investigate conditions not readily possible with ground or flight test facilities. The solutions presented in the 
following were computed with the Reynolds Averaged Navier-Stokes equations using the Spalart-Allmaras 
turbulence model. The limitations of RANS for post-stall flows are well known. One objective of these 
computations is to assess the applicability of the RANS methodology in modeling the aggregate dynamic behavior 
of transport aircraft into and through stall. All solutions to follow, both static and dynamic, were advanced with 
second order physical time stepping. 

Figure 8 provides a comparison of wind tunnel and CFD rolling-moment coefficient results for a range of angle 
of attack and frequency. These results correspond to GTM roll-forced sinusoidal oscillations with amplitude 20°. 
The oscillations were initiated from converged time-accurate static solutions at each angle of attack, and rotated for 
two complete oscillations. The dynamic solutions were essentially converged after one-fourth of a cycle. The/=0.92 
Hz solution was advanced with 360 time steps per roll cycle, and th e/ = 0.12 and 0.06 Hz cases with 1800 time 
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steps per cycle. The red lines on Fig. 8 indicate the second roll cycle of the CFD simulation and black corresponds to 
the wind- tunnel mean-cycle measurements. Light gray lines show the wind-tunnel measurements before averaging 
the cycles of oscillation. These measurements highlight the substantial measurement noise and other responses 
outside the desired bandwidth at high angles of attack where separated flows dominate. Overall, a reasonable match 
in response shape is presented although at a 0 = 12° some differences occur reflecting the sensitivity of modeling this 
difficult unsteady flight condition. 

To further explain the modeling sensitivity near the unstable and unsteady condition at a 0 = 12°, Fig. 9 presents 
streamlines over the wing, that were computed to show how the flow is changing around this flight condition. These 
static solutions were advanced in time at a characteristic time step of 0.02, i.e. a free stream particle will travel the 
length of the mean aerodynamic chord in 50 time steps. In addition to streamlines, lift and pitching moment static 
curves are presented below the wings along with the corresponding wind-tunnel data. Figure 9 shows a sequence of 
streamlines over the wing for four flight conditions: Fig. 9a, a = 4°, showing normal flight conditions; Fig. 9b, a = 
11.75°, showing conditions just before a break in the lift curve; Fig. 9c, a - 12°, showing stalled conditions; and 
Fig. 9d, a - 12.5°, showing post stall-break conditions. The figure shows that when the abrupt change in lift occurs 
at a - 12°, streamlines also make a rapid change and indicate a stalled wing tip. Also in this region, the lift curve 
slope has become flat on either side of the break in lift. These adverse conditions support the unstable roll damping 
found in forced oscillation testing (Fig. 4) and is confirmed by abrupt roll departures demonstrated in GTM subscale 
flight tests near this angle of attack. CFD results point to a more abrupt change in lift and pitching moment than 
found in the wind-tunnel data. 

CFD simulation was also used to generate more detailed forced-oscillation data over a limited range of a 0 in key 
regions of interest. Two regions in a 0 were computed that covered 10° < a 0 < 15° and 24° < a 0 < 30° with fine 1° 
increments. Figure 10 presents harmonic analysis for CFD simulated roll forced-oscillation data. To reduce clutter, 
only high frequency results from the wind-tunnel analysis are overlaid for comparison. In general, the CFD analysis 
produced larger magnitudes for the coefficients but identified similar unsteady regions. The plots show that CFD 
results imply a shift of the largest roll instability toward a 0 = 14° rather than at a 0 = 12° as indicated by the wind- 
tunnel results. CFD results indicate unsteady behavior over a wider range of 11° < a 0 < 15°, compared to the wind- 
tunnel result implying just at a 0 = 12°. In the lower a 0 range R 2 analysis is similar to the wind-tunnel results. In that 
range both wind-tunnel and CFD capture the nonlinearity at high frequency. At the upper a 0 range, the CFD results 
present higher values of R 2 reflecting the lower contribution of noise and other higher order effects in CFD 
simulation and confirming the use of linear unsteady models. 

D. Unsteady Model Parameter Estimation 

Parameter estimation was accomplished using a variation on an output-error (OE) method 7,8,20 in the time 
domain. In the OE approach the model equations are given by Eqs. (9) and (10). Measured inputs, P , and outputs, 
C a , for each of the frequencies tested, are combined to ensure all the frequency content is included in the time 
domain estimation process. General application of the output-error method for aircraft is explained and relevant 
software is provided in Ref. 2 1 . 

Results from application of OE to GTM roll forced-oscillation data is shown in Fig. 11. This figure shows the 
estimated parameters and their 2a bounds for a large range of angle of attack, for 0° < a 0 < 70°. Uncertainty bounds 
are generally very small and are in statistical agreement. This implies good information content or that adequate 
sensitivity was obtained for parameter estimation using the wind-tunnel data. The top plot shows steady-flow 
damping. Consistent with previous results, this term is neutral to unstable at a 0 - 12° and a 0 > 50°. The next two 
plots, in Fig. 11, define the unsteady term. The middle plot shows model parameter, a , that acts as a gain on the 
unsteady term and is, in general, relatively small except at a 0 = 12°. The bottom plot presents a non-dimensional 
time constant defining lag in the aerodynamic response. This plot shows three regions where unsteady behavior is 
present and the relative size of the lags in those regions. The smallest lag or fastest response is at a 0 = 12°. 

As an example of model fidelity, at a 0 = 30°, response prediction is presented in Fig. 12. This example shows 
measured and unsteady model responses for both low and high frequency cases of the roll moment against sideslip. 
Response measurements plotted against ^ are commonly used and have been presented in this study. However for 
this graphic, using p highlights the adequacy of the linear unsteady model even with relatively large displacements 
of this key state variable. Equation (6) defines the relationship between $ and p. The estimated linear static curve is 
included to provide a reference line. Since linear models are used, the computed responses are regular ellipses. In 
general the models represent the measured data very well. One notable feature for the low frequency case is the very 
poor signal-to-noise for this test point. At the higher frequency, signal-to-noise is improved but still not very good. 
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The measured data for the high frequency case does reflect some nonlinearity as indicated by the slightly irregular 
elliptical shape of the data; however, the linear model is still an adequate representation in this case. Although lag in 
the response caused by the unsteady flow is not readily visible in this type of plot, some evidence is presented by the 
rotation of the ellipse relative to the static line. For a linear model without a lag in the response, the long axis of the 
ellipse will align with the static line. 


VI. Concluding Remarks 

This paper presents a study of roll-axis aerodynamics for the NASA GTM subscale vehicle. A method for 
harmonic analysis is presented and was performed on both experimental and computational force-oscillation data to 
identify where unsteady models are required and to characterize the linear or nonlinear nature of the models. 
Harmonic analysis was also applied to estimate cross derivatives and to determine the impact of amplitude effects. A 
best practice was suggested for harmonic analysis and R 2 analysis by appropriately filtering measurement data. A 
more detailed investigation of roll-axis instabilities was performed. The approach for investigating difficult 
aerodynamic flight regimes in aircraft suggested using both experimental and analytical technologies to provide 
greater insights into the aerodynamic phenomenon studied. Wind-tunnel dynamic data and simulated computational 
data showing reasonable comparisons support this approach. Computational methods provided aerodynamic details 
(dynamic force and moment data and streamlines) over a fine angle -of-attack grid near unsteady and unstable 
conditions. Computational data presented suggest the presence of roll instability (in agreement with the experimental 
analysis) and a lift curve stall break at 12° angle of attack. The abrupt break in lift was not observed in the wind- 
tunnel results; however, flight testing does support the presence of an abrupt roll instability at this condition. Using 
the output-error estimation method, linear unsteady models were identified for the roll moment that can support a 
global nonlinear representation. An example of estimated model fidelity was demonstrated at a high angle of attack, 
a 0 = 30°, where unsteady aerodynamic behavior was found. The unsteady model represented the measured 
responses well. 

Application of the unsteady modeling methodology to static and dynamic wind-tunnel data allows simulations of 
transport aircraft responses in adverse conditions to be developed and, in turn, facilitate development of control 
designs that provide stability and control in difficult flight regimes. 
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Figure 1. Model geometry for NASA 
experimental subscale aircraft. 
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Figure 2a. Harmonic analysis for rolling moment 
coefficient, roll oscillations, <j> A = 10°. 
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Figure 2b. Harmonic analysis for rolling moment 
coefficient, roll oscillations, (j) A = 20°. 
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Figure 2c. Harmonic analysis for rolling moment 
coefficient, roll oscillations, ^ 4 = 30°. 
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Figure 3a. Harmonic analysis for yawing 
moment coefficient, roll oscillations, </> A = 20°. 
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Figure 3b. Harmonic analysis for side 
force coefficient, roll oscillations, (j) A = 20°. 
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Figure 4a. Harmonic analysis for rolling moment 
coefficient, roll oscillations, (j) A = 20° (4 Hz filter). 




Figure 4b. Mean-cycle and first order harmonic 
model responses for rolling moment coefficient, 
roll oscillations, a 0 = 30°, (/> A = 20°. 
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Figure 5a. Harmonic analysis for side force 
coefficient, roll oscillations, (j> A = 20° (10 Hz filter). 
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Figure 5b. Harmonic analysis for side force 
coefficient, roll oscillations, (f> A - 20° (2 Hz filter). 
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Figure 6a. Harmonic analysis for rolling moment 
coefficient, roll oscillations, (j> A - 20° (10 Hz filter). 
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Figure 6b. Harmonic analysis for rolling moment 
coefficient, roll oscillations, (j) A - 20° (2 Hz filter). 
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Figure 8. Rolling moment forced oscillations from CFD and wind tunnel (unfiltered and mean-cycle) 
data, with amplitude (/> A = 20°, for GTM over range of angle of attack and frequency. 
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Figure 9a. CFD calculation of streamlines 
for GTM at a = 4°. 


Figure 9b. CFD calculation of streamlines 
for GTM at a= 11.75°. 
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Figure 9c. CFD calculation of streamlines 


Figure 9d. CFD calculation of streamlines 


for GTM at a =12°. 


for GTM at a =12.5°. 
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Figure 11. Unsteady model parameters with 
2a uncertainty bounds for rolling moment 
coefficient, </> A =20°. 


Figure 10. Harmonic analysis for CFD simulated 
rolling moment coefficient, roll oscillations, </> A = 20° 
(4 Hz filter applied to wind tunnel data). 
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Figure 12. Unfiltered measurements and unsteady model responses for rolling moment 
coefficient, at two oscillation frequencies, a= 30°, (j) A = 20°. 
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